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Abstract 

We consider the heat equation in the presence of compactly supported 
magnetic field in the plane. We show that the magnetic field leads to 
an improvement of the decay rate of the heat semigroup by a polynomial 
factor with power proportional to the distance of the total magnetic flux 
to the discrete set of flux quanta. 

The proof employs Hardy-type inequalities due to Laptev and Weidl 
for the two-dimensional magnetic Schrodinger operator and the method of 
self-similar variables and weighted Sobolev spaces for the heat equation. 
A careful analysis of the asymptotic behaviour of the heat equation in the 
similarity variables shows that the magnetic field asymptotically degener- 
ates to an Aharonov-Bohm magnetic field with the same total magnetic 
flux, which leads asymptotically to the gain on the polynomial decay rate 
in the original physical variables. 

Since no assumptions about the symmetry of the magnetic field are 
made in the present work, it confirms that the recent results of Kovaf fk [TJ] 
about large-time asymptotics of the heat kernel of magnetic Schrodinger 
operators with radially symmetric field hold in greater generality. 
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1 Introduction 



The principal objective of this paper is to show that the solutions of the heat 
equation 

u t + (-iV x -A{x)) 2 u = 0, (1) 

in space-time variables (x, t) £ R 2 x (0, oo), admit a faster large-time decay due 
to the presence of a magnetic field B — not A that is supposed to vanish at 
infinity. The characteristic property under which we prove the result is that the 
total magnetic flux 

$ s :=±- [ B(x)dx (2) 

does not belong to the discrete set of flux quanta, which coincides with integers 
for our choice of physical constants. As a measure of the additional decay of 
solutions, we consider the (polynomial) decay rate 

1B :=sup{ 7 | 3C 7 >0, Vt >0, \\e~ { ^ V ~ A)2t \\ L2{K) ^ L2 <C 7 (l+t)-^ . 

(3) 

Here we use the shorthands L 2 := L 2 (R 2 ) and L 2 (K) := L 2 (R 2 , K(x) dx) with 
the Gaussian weight 

K(x) := el x l 2 / 4 . (4) 
Our main result reads as follows: 
Theorem 1. Let B £ CSfR 2 ). We have 



IB 



= 1/2 ifB = Q, 

> (1 + ft)/ 2 otherwise, 



where ft := dist($s, 



Here the result for B = follows easily by the explicit knowledge of the 
heat kernel for the Laplacian in R 2 . The essential content of Theorem [1] is 
that solutions to ((T|) when the total flux $s is not an integer gain a further 
polynomial decay rate /3/2. The proof of this statement is more involved and 
constitutes the main body of the paper. It is based on the method of self-similar 
solutions recently developed for twisted-waveguide systems by Zuazua and the 
present author in [161 115] . 

The result has a rather natural interpretation in terms of diamagnetism that 
we briefly explain now. 



1.1 Relation to Hardy inequalities 

Diamagnetism in quantum mechanics is reflected in repulsive properties of the 
Hamiltonian 

H B := HV - A) 2 (5) 

describing a spin- less particle in the external magnetic field B. This is most 
easily seen for the homogeneous field B{x) = Bq for which the ground-state 
energy is raised by the field magnitude, i.e. inf cf(Hb q ) = \Bq\, and any such non- 
zero field therefore leads to a Poincare-type inequality Hb > \B \. The latter 
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remains true for magnetic fields satisfying ±B(x) > |-Bo|- This clearly implies an 
exponential decay rate \\e~ HBt \\ L -2^ L 2 < e~l B °l* in view of the spectral mapping 
theorem. 

For magnetic fields vanishing at infinity considered in this paper, the repul- 
sive effect is more subtle because the spectrum starts by zero: ct(Hb) = [0, oo), 
regardless of the presence of the magnetic field. Hence, \\e~ HBt \\L 2 ^L 2 — 1- 
However, Laptev and Weidl proved in |17] an important Hardy-type inequality 

in the sense of quadratic forms, where c = c(A) is a positive constant whenever 
<&b ^ Z (for generalizations, see [2TJ 121 HH] ) ■ On the other hand, because of the 
criticality of the free Hamiltonian in R 2 , ([6]) cannot hold for A = 0. It is natural 
to expect - and supported by the general conjectures raised in [THl H2] - that 
the additional repulsiveness caused by the Hardy-type inequality will lead to an 
improvement of the asymptotic decay rate of solutions to ((T|) , provided that the 
initial data are restricted to a subspace of L 2 . In fact, ^ plays a central role 
in our proof of Theorem [TJ 

For dimensions greater or equal to three, the inequality © is just a di- 
rect consequence of the diamagnetic inequality (|35l) and the classical Hardy 
inequality for the free Hamiltonian, so that ([5]) holds with a positive constant 
independent of A in the higher dimensions. That is why the most interesting 
situation is the two-dimensional case, to which we restrict in this paper. 



1.2 Relation to previous results of Kovafik |14| 

Before elaborating on the details of the proof of Theorem [TJ let us make a few 
comments on the novelty of the results. 

The present paper is partially motivated by a recent study of Kovafik [14] 
about large-time properties of the heat kernel of two-dimensional Schrodinger 
(and Pauli) operators with radially symmetric magnetic field. Among other 
things, the author of [14] shows that the magnetic heat kernel has a faster decay 
for large times than the heat kernel of the free Hamiltonian. In some cases, he 
even shows that the additional decay rate is exactly that of Theorem [TJ His 
results are more precise in the sense that he gets pointwise (sometimes two- 
sided) estimates and cover some cases with $s € Z. 

On the other hand, the method of [M] is based on partial- wave decomposition 
and is therefore restricted to radially symmetric fields. The method of the 
present paper, on the contrary, is more general and it is actually relying only on 
the existence of the magnetic Hardy- type inequality ([Fj]) . In this way, we have 
been able to treat the general magnetic field without any symmetry restriction. 

Finally, Kovafik calculated in [14] the exact large-time asymptotics for the 
Aharonov-Bohm field. These asymptotics yield the same decay rates as for 
regular magnetic fields, suggesting a somewhat universal nature of the decay 
rates. The present paper provides an insight into the phenomenon: all the 
magnetic fields effectively behave for large times as an Aharonov-Bohm field. 
This statement will be made precise in the following subsection. 
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1.3 The idea of the proof 



As already mentioned above, the principal idea behind the main result of The- 
orem [TJ i.e. the better decay rate due to the presence of magnetic field, is the 
positivity of the right hand side of (|6j) . In fact, Hardy- type inequalities have 
already been used as an essential tool to study the asymptotic behaviour of the 
heat equation in other situations [5J [TH1 [H] . However, it should be stressed 
that Theorem [1] does not follow as a direct consequence of (JSJ) by some energy 
estimates but that important and further technical developments that we ex- 
plain now are needed. Nevertheless, overall, the main result of the paper again 
confirms that the Hardy-type inequalities end up enhancing the decay rate of 
solutions. 

Let us now informally describe our proof that there is the extra decay rate 
if the magnetic field is present. 

I (changing the space-time). The main ingredient in the proof is the method 
of self-similar solutions developed in the whole Euclidean space by Escobedo and 
Kavian [8] . We rather follow the approach of [16l [15] where the technique was 
adapted to twisted-waveguide systems which exhibit certain similarities with 
the present problem. Writing 

u(y,s)^e s /Me s/2 y,e s ~l), (7) 

the equation (TTJ) is transformed to 

u s -\y-V y u+ (-iVj, - A s {y)) 2 u-\u = (8) 

in "self-similarity" variables (y,s) £l 2 x (0, oo), where 

A s (y)--=e s/2 A(e s/2 y)- (9) 

When evolution is posed in that context, y plays the role of new space variable 
and s is the new time. However, note that now we deal with a non-autonomous 
system. 

II (choosing rapidly decaying initial data). We reconsider © in the 
weighted space L 2 (K). It has the advantage that the associated (s-dependent) 
generator has compact resolvent in this space and therefore purely discrete spec- 
trum. The latter is particularly useful for energy estimates. 

III (removing the weight). Just for personal reasons, we prefer to work 
in the usual I/ 2 -space, without the weight K. To achieve this, we make the 
substitution v — K x l 2 u and transform © on L 2 (K) to the equivalent problem 

v s +(-iV y -A s (y)) 2 v+-^-v = (10) 

on L 2 . We observe the presence of the harmonic-oscillator potential that makes 
the spectrum of the generator to be discrete. Furthermore, recalling the spec- 
trum of the quantum-harmonic-oscillator Hamiltonian (see any textbook on 
quantum mechanics, e.g., [131 Sec. 2.3]) and the relationship between s and t 
from (J7|), it is clear now why we get the 1/2 for the decay rate in Theorem [T] if 
A = 0. It remains to analyse the effect of the presence of the magnetic field. 

IV (taking s — > oo). We look at the asymptotic behaviour of (|10p as the 
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self-similar time s tends to infinity. The scaling coming from the self-similarity 
transformation is such that 

B s :=rotA s >2n$ B 8=:B 00 (11) 

in the sense of distribution, where 8 is the Dirac delta function and <&b is total 
flux (0) . The right hand side of (fTTj) is the Aharonov-Bohm field to which there 
corresponds an explicit vector potential having a singularity at 0. Taking 
this limit into account, it is expectable that (|10[) will be approximated for large 
times by the Aharonov-Bohm problem 

w s + (-iV v -A QC (y)fw + -^jrW = 0. (12) 

This evolution equation is explicitly solvable in terms of special functions. In 
particular, it is clear that the solutions decay exponentially with the decay 
rate given by the lowest eigenvalue of the (time-independent) generator which 
explicitly reads 

(l+/3)/2. (13) 

It remains to observe that the exponential decay rate transfers to a polynomial 
one after returning to the original time t and that the number (|13p is just that 
of the bound of Theorem [1] in the presence of magnetic field. 

Two comments are in order. First, although it is very natural to expect, 
we do not establish any theorem that solutions of (|10j) can be approximated 
by those of (|T2|) as s — > oo. We only show a strong-resolvent convergence for 
operators related to their generators (Proposition^). This is, however, sufficient 
to prove Theorem [T] with help of energy estimates. Proposition 0] is probably 
the most significant auxiliary result of the paper and we believe it is interesting 
in its own right. Second, in the proof of Proposition 0] we heavily rely on the 
existence of the Hardy- type inequality ((6]). 

1.4 The content of the paper 

The organization of this paper is as follows. 

In the following Section [5] we give a precise definition of the magnetic Schro- 
dinger operator ([5]) and the associated semigroup. We also briefly discuss there 
the case of zero magnetic field, obtaining, inter alia, the first statement of 
Theorem Q] ( cf Proposition [1] and Corollary [1} . 

The main body of the paper is represented by Section [3] where we develop 
the method of self-similar solutions to get the improved decay rate of Theo- 
rem [T] as described above. Here we point out Section 13.41 where a convergence 
of operators associated with the regular and Aharonov-Bohm fields of (jll|) is 
proved (cf Proposition 01 . Finally, in Section |3~51 we establish an alternative 
version of Theorem [TJ which provides a global (in time) upper bound to the heat 
semigroup (cf Theorem [2]). 

The paper is concluded in Section 0] by referring to some open problems. 

2 Preliminaries 

Let the local magnetic field B be represented by a real-valued continuous func- 
tion of compact support in K 2 . It is well known that the corresponding magnetic 
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vector potential A is not uniquely determined by the formula 



B = rot A = dxA 2 - d 2 A 1 . (14) 

This freedom, reflecting the gauge invariance of the physical theory, enables us 
to work with the particular choice 



A{x) := (— X2,X\) / B[txx,tX2)t dr (15) 
Jo 

without loss of any generality. It corresponds to the so-called multipolar gauge 
in molecular electrodynamics, and it is also known as the transversal gauge 
because x ■ A(x) = 0. For smooth B it is straightforward to verify the validity 
of HU). If B is merely continuous, the identity (fT4"]) can be checked in the sense 
of distributions by approximating B by a sequence of mollifiers. 

In polar coordinates (r, 9) G [0, oo) x [0, 2ir), related to the Cartesian system 
by x — (r cos 9, r sin 9), the magnetic potential (|15[) acquires a more transparent 
form 

A(x) = (-sin0, cos 9) (16) 

r 

with the notation 



a(r, 9) 



I B(t cos 9,t sin 9) t dr 
Jo 



Since B is supposed to be compactly supported, the function a is bounded 
and A is vanishing at infinity. Moreover, the limit 

aoo (9) := lim a(r,9) (17) 

is well defined for every 9 <E [0, 2tt). 
The quantity 



1 f 27T 
$(r) := — / a(r,9)d9 
27r Jo 



has the physical meaning of the magnetic flux through the disc of radius r 
centred at the origin of M 2 . Since the field B is zero outside such a disc of 
sufficiently large radius, $(r) = <&b for all r large enough, where $s is the total 
magnetic flux introduced in @. As pointed out in the introduction, the set of 
integers correspond to flux quanta for our choice of physical constants (cf |131 
Sec. 10.2.4]). 

The magnetic Hamiltonian (JS|) is defined as the self-adjoint operator associ- 
ated on L 2 :— L 2 (K 2 ) with the closed quadratic form 

h B [$\ := \\Vij~iA^\\ 2 , ^S(ft B ):=^ 2 , (18) 

where || • || denotes the norm in L 2 and W 1 ' 2 :— M /1,2 (R 2 ) is the usual Sobolcv 
space. Using the transverse gauge and polar coordinates, we can write 

h B m = \\d r ^\\ 2 + \\r-\deil> - iai>)\\ 2 ■ (19) 

Remark 1. By the gauge invariance mentioned above, it makes sense to dis- 
tinguish the Hamiltonian by the subscript B rather than A. Indeed, all the 
Hamiltonians corresponding to different gauges are unitarily equivalent (c/ |13l 
Prob. 4.531). 
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Since Hb — i?o is a relatively form compact perturbation [HO Prob. XIII. 39] 
of the free Hamiltonian Hq, we have 

a(H B ) = [0,oo). (20) 

As usual, we consider the weak formulation of the parabolic problem ^ 
subjected to the initial datum uo eL 2 . We say a Hilbert space- valued function 
u e Lf oc ((0, oo); W 1 ' 2 ), with the weak derivative v! € L 2 OC ((0, oo); TU" 1 ' 2 ), is a 
(global) solution of (HJ provided that 

<u,u'(t)> + /i B ( 1 ;,u(t)) =0 (21) 

for each v £ W 1 ' 2 and a.e. i e [0, oo), and zt(0) = tio- Here /is(-, •) denotes the 
sesquilinear form associated with (IT%)) and (•, •) stands for the pairing of W 1 ' 2 
and its dual W~ 1,2 . With an abuse of notation, we denote by the same symbol u 
both the function on M 2 x (0, oo) and the mapping (0, oo) — > W 1 ' 2 . 

Standard semigroup theory implies that there indeed exists a unique solution 
of (HH) that belongs to C°([0, oo), L 2 ) . More precisely, the solution is given 
by u(t) = e~ HBt Uo, where e~ Hmt is the semigroup associated with Hb- It is 
important to stress that the semigroup is not positivity-preserving, which in 
particular means that we have to work with complex functional spaces when 
studying (121]) . 

The spectral mapping theorem together with (|20l) yields \\e~ HBt \\ L 2^ L 2 = 1 
for each time t > 0. To observe an additional large-time decay of the solu- 
tions u(t), we have to restrict the class of initial data uo to a subspace of L 2 . In 
this paper we consider the weighted space L 2 (K) := L 2 (M. 2 , K(x) dx) with the 
weight given by ((4]). The norm in L 2 (K) will be denoted by |j • \\k- 

The case of the semigroup of the free Hamiltonian Hq is well understood 
because of the explicit knowledge of its heat kernel: 

Using this expression, it is straightforward to establish the following bounds: 



Proposition 1. There exists a constant C such that for every t > 1. 

C-H-" 2 <\\e- H ^\\ L , {K) ^ <Ct-" 2 . 
Proof. Using the Schwarz inequality, we get 

IK^oll 2 < \\u f K I lMX K *'f 2 dxd X > 

for every uo G L 2 (K), where the integral can be evaluated explicitly and yields 
the upper bound. To get the lower bound, it is enough to find one particular 
u 6 L 2 (K) such that ||e~ Hot u || > C^ 1 ^ 1 / 2 \\u \\ K . It is just a matter of 
explicit computations to check that it is the case for the choice uq = with 
any j3 > 1/2. □ 

As a consequence of this proposition, we get: 

Corollary 1. We have 70 = 1/2. 

This establishes the trivial part of Theorcm[TJ On the other hand, no explicit 
formula for the heat kernel is available if B ^ and a more advanced technique 
is needed to study the decay rate 73. 
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3 The self-similarity transformation 



Our method to study the asymptotic behaviour of the heat equation |T]) in the 
presence of magnetic held is to adapt the technique of self-similar solutions used 
in the case of the heat equation in the whole Euclidean space by Escobedo and 
Kavian [5] to the present problem. We closely follow the approach of the recent 
papers [161 115) , where the technique is applied to twisted waveguides in three 
and two dimensions, respectively. 

3.1 An equivalent time-dependent problem 

We consider a unitary transformation U on L 2 which associates to every solution 
u £ L 2 OC ((0, oo), dt; L 2 ) of (j2Tj) a function u = Uu in a new s-time weighted 
space L 2 oc ((Q,oo),e s ds;L 2 ) via ©. The inverse change of variables is given by 

u(x, t) = (t + I)- 1 ' 2 u((t + l)- 1/2 x, log(t + 1)) . (22) 

Note that the original space-time variables (x, t) are related to the "self-similar" 
space-time variables (y, s) via the relations 

(x,t) = {e s ' 2 y,e s - 1) , (y, s) = ((t + iy 1/2 x, log(< + 1)) . (23) 

Hereafter we consistently use the notation for respective variables to distinguish 
the two space-times. 

It is easy to check that this change of variables leads to the evolution equa- 
tion ([8]), subject to the same initial condition as u in ([1]). More precisely, the 
weak formulation l|21|) transfers to 

(v,u'(s)-±y V v u{s)) + a s (v, u{s)) = (24) 

for each v £ W 1 ' 2 and a.e. s £ [0, oo), with w(0) = uq := Uuq = u . Here a s (-, •) 
denotes the sesquilinear form associated with 

a.[u] := ||Vtt - iA s u\\ 2 - i ||ii|| 2 , u G := W 1 ' 2 , 

where A s has been introduced in (]9]). 

Remark 2. Note that |5]) (or is a parabolic equation with s-time-dependent 
coefficients. The same occurs and has been previously analysed in the twisted 
waveguides [16l [15] and also for a convection-diffusion equation in the whole 
space but with a variable diffusion coefficient [7] . A careful analysis of the 
behaviour of the underlying elliptic operators as s tends to infinity leads to a 
sharp decay rate for its solutions. 

Since U acts as a unitary transformation on L 2 , it preserves the space norm 
of solutions of ([T]) and dH), i.e., 

Kt)|| = ||fi(*)||. (25) 

This means that we can analyse the asymptotic time behaviour of the former 
by studying the latter. 

However, the natural space to study the evolution © is not L 2 but rather 
the weighted space L 2 (K). Following the approach of 16 based on a theorem 
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of J. L. Lions [2J Thm. X.9] about weak solutions of parabolic equations with 
time-dependent coefficients, it can be shown that ([5]) is well posed in the scale 
of Hilbcrt spaces 

W 1>2 (K) C L 2 {K) C W 1>2 (K)*, (26) 

with 

W h2 (K) := {u e L 2 (K) | Vu G L 2 {K)} . 

More precisely, choosing v := Kv for the test function in ([24|. where v £ 
C(j°(rio) is arbitrary, we can formally cast (|2"4")l into the form 

(u,{^s))^ + a s («,i/(s)) =0. (27) 

Here (•, denotes the pairing of W 1,2 {K) and W 1,2 (K)* , and 

a a (u,it) := (Vu - i^4 s -y, Vu - iVU/)^. - - (yu, L4yit) K - - , (28) 

with (-,-)^f denoting the inner product in L 2 (K). Note that a s is not a sym- 
metric form. 

By "formally" we mean that the formulae are meaningless in general, be- 
cause the solution u(s) and its derivative w'(s) may not belong to W l ' 2 (K) 
and W 1,2 (K)* , respectively. The justification of © being well posed in the 
scale (|26| consists basically in checking the boundedness and a coercivity of the 
form a s defined on D(a s ) := W 1 ' 2 ^). It is straightforward by noticing that A s 
is bounded for each s fixed. We refer to [TH1 Sec. 5.4] for more details. 

3.2 Reduction to a spectral problem 

Choosing v := ii(s) in (1271) and combining the obtained equation with its con- 
jugate version, we arrive at the identity 

~\\Hs)f K = -i s [u(s)}, (29) 

where l s [u] := 3?{a s [{(]}, u 6 S)(Z S ) := D(a s ) = W 1 ' 2 ^) (independent of s). 
Observing that the real part of the non-symmetric term in (1281) vanishes, we get 

l„[u] := \\Vu - iA s u\\ K - - \\u\\ 2 K . 

It remains to analyse the coercivity of l s . 

More precisely, as usual for energy estimates, we replace the right hand side 
of (|2"9")l by the spectral bound, valid for each fixed s 6 [0, oo), 

Vu l s [u]>X(s)\\u\\ 2 K , (30) 

where A(s) denotes the lowest point in the spectrum of the self-adjoint opera- 
tor L s associated on L 2 (K) with l s . Then (|2T)1) together with ([3^| implies the 
exponential bound 

V S £[0,oo), \\u(s)\\ K < \\u \\ K e -/o^)rfr_ (31) 

Finally, recall that the exponential bound in s transfers to a polynomial bound 
in the original time t, cf (|23p . In this way, the problem is reduced to a spectral 
analysis of the family of operators {L s } s >o. 
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3.3 Study of the spectral problem 

In order to investigate the operator L s in L 2 (K), we first map it into a unitarily 
equivalent operator L s :— ULJJ^ 1 on L 2 via the unitary transform 

Uu := K 1/2 u. 

By definition, L s is the self-adjoint operator associated on L 2 with the quadratic 
form l s [v] :— l s [U~ 1 v], v e 25(/ s ) := UQ(l s ). A straightforward calculation 
yields 

l,[v] = ||Vw - iA s v\\ 2 + ^ IMI 2 , v G 3)(i.) = W X ' 2 (R 2 ) n L 2 (R 2 , |2/| 2 dy) . 

(32) 

In particular, 2)(Z S ) is independent of s. In the distributional sense, we can write 

L s = {-iV-A s f + ±-\y\ 2 . (33) 
16 

For the trivial magnetic field B = 0, the operator L s coincides with the 
Hamiltonian of the quantum harmonic oscillator 



1 

16 



, HO := -A + i- \y\ 2 in L 2 (R 2 ) (34) 



(i.e. the Friedrichs extension of this operator initially defined on C^°(1R 2 )). In 
any case, the form domain of Lho coincides with 23 (l s ). Recall the well known 
fact that Lho has compact resolvent (see, e.g., [HI Thm. XIII. 67]) and thus 
purely discrete spectrum. We conclude that, for any field B, L s has compact 
resolvent, too. In particular, A(s) represents the lowest eigenvalue of L s . 
By the diamagnetic inequality [TT] Thm. 2.1.1] 

| (— iV — A s )v\ > \V\v\\ (35) 

valid for every v G L 2 almost everywhere and by the minimax principle, it 
follows that A(s) is bounded from below by the lowest eigenvalue of Lho- Using 
the explicit knowledge of the spectrum of Lho (see, e.g., [131 Sec. 2.3]), we get 

VsG[0,oo), A(s)>1/2. (36) 

Obviously, there is an equality if B = 0. Moreover, the bound becomes sharp 
in the limit s — > oo for magnetic fields with the total flux being an integer. 

Proposition 2. Let <f> B G Z. Then 

A(oo) := lim \{s) = 1/2. 

s— >oo 

Proof. In view of f|36[) . it is enough to establish an upper bound to A(s) that 
goes to 1/2 as s — > oo. We obtain it by constructing a suitable trial function in 
the variational characterization of the eigenvalue. 

For every n G N larger than 1, let us define the function 

{1 if r > 1/n , 

if r G [1/n 2 , 1/n] , 
logn 
otherwise . 
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We have \r/ n \ < 1 for all n > 2, r/ n (r) converges as n — >• oo to 1 for every 
r G (0, oo) and 

WJh«0,oo),Tdr) = ( lo g«) _1 ► ' (37) 



Hence, the functions (r, 9) H> r) n (r) in polar coordinates in R 2 represent a conve- 
nient approximation of the constant function 1 in W lo ' c (M 2 ) when a singularity 
in the origin has to be avoided. 

In polar coordinates in R 2 , we then introduce 



where (r, 6) i— > := e~ r / 8 is an eigenfunction of Lho corresponding to its 
lowest eigenvalue 1/2 (cf Proposition [3]) and is defined in (fl7|) . Since ^„ 
vanishes in the vicinity of the origin and decays sufficiently fast at the infinity, 
we clearly have ip n G r D(l s ) for all n > 2 and s > 0. The dominated convergence 
theorem implies 



llV'nll > 2tt / (f{r) z r dr — An . 

n->oo J 

At the same time, a straightforward calculation yields 

/>oo />2tt 

UV>»] - lll^f = / / ip(r) 2 r, n (r) 2 \ aoo (e)-a(e s / 2 r,e)\ 2 de '- 
Jo Jo 

/•OO 

+ 2n ip(r) 2 rj' n (r) 2 r dr . 



Since a(e s / 2 r, 0) converges as s — > oo to for every (r, 9) G (0, oo) x [0, 2-7r) 
and a is a bounded function, the double integral tends to zero as s — > oo for every 
n > 2 fixed due to the dominated convergence theorem. The one-dimensional 
integral is independent of s and it goes to zero as n — > oo due to (|37[) and 

M<i. 

By (j3l)|) and by the Rayleigh-Ritz variational formula for A(s), we have the 
bounds 

2 WnW 

for every s > 0. Fixing n > 2 and taking the limit s — > oo, we therefore get 

< A(oo) - - < . 

Since the right hand side provides an upper bound which can be made arbitrarily 
small by choosing n sufficiently large, we conclude with the statement of the 
proposition. □ 

Proposition [5] explains why the estimate (|3ip does not lead to any improved 
polynomial decay rate if $s G Z. To get the improved decay rate for <&b & 
Z, we need to show that a better inequality than (|3l)f is available, at least 
asymptotically as s — > oo. 



11 



3.4 The asymptotic behaviour of the spectrum 

This subsection is devoted to a study of the asymptotic behaviour of the oper- 
ator L s as s — > oo. Since the limit directly leads to the improved decay rate, it 
is probably the most essential part of the paper. 

We recall the expression (|16l) for the magnetic potential A in the transverse 
gauge and conventionally expressed in polar coordinates. It is clear that the 
scaling © acts in such a way that 

A s > locally uniformly in R 2 \ {0} , (38) 

s— ► oo 

where (recall fT7])) 

A^iy) := (- S in0,cos0) (39) 

r 

is well defined off the origin. Here, to save symbols, we use the same notation 
(r, 9) G [0, 00) x [0, 2ir) for polar coordinates in the "self-similar" y-plane as well. 
Contrary to A s , which is a bounded function for every s G [0, 00), A^o possesses a 
singularity at the origin. As already revealed in the introduction, A^ represents 
the magnetic potential in the transverse gauge for the Aharonov-Bohm field 
defined in (fTTj) . Indeed, rot Aoo = 2tt&b S in the sense of distributions. 

3.4.1 The Aharonov-Bohm Hamiltonian 

In view of ((38)) . it is reasonable to expect that L s converges in a suitable sense 
as s — > 00 to the operator 

Loo := (-*V - A^) 2 + -L \y\ 2 . (40) 
Id 

However, one has to be careful when giving a meaning to Loo because of the 
singularity of Aoo- Indeed, already for the operator (— iV — Aoo) 2 with a non- 
trivial rotationally-symmctric Aharonov-Bohm field, i. e. Qk, (9) — const ^ Z, it 
is well known that the symmetric operator defined by this differential expression 
on the minimal domain C ( J o (R 2 \{0}) is not essentially self-adjoint and that there 
are infinitely many self-adjoint extensions determined by imposing boundary 
conditions at the origin (cf [18] for a modern treatment in terms of boundary 
triplets) . 

It turns out that the "limit" of L s as s — > 00 corresponds to the Friedrichs 
extension of Loo, *-e-, we understand (|40[) as the self-adjoint operator associated 
on L 2 with the quadratic form 

U«] := || Vu - iA^vf + i ||y V || 2 , v G D^) := C °°(IR 2 \ {O})"'"' 00 . 

(41) 

Here the norm || ■ 11;^ with respect to which the closure is taken is defined by 
IM| 2 : = loo[v] + |M| 2 - We call Loo simply the Aharonov-Bohm Hamiltonian, 
although it contains the additional quantum-harmonic-oscillator potential term 
with respect to the classical Aharonov-Bohm Hamiltonian. 

Now we specify the domain of the form loo- Observing the structure of loo in 
polar coordinates (cf (|19l) ) and being inspired by [17] . we consider the ordinary 
differential self-adjoint operator K on the Hilbert space L 2 ((0, 27r)) defined by 

K<t>:=i(f>' + a 00 </>, 0eD(X) := {0 G W 1 ' 2 ^, 2tt)) | 0(0) = 0(2vr)} . 
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The spectral problem for K is explicitly solvable. In particular, we have <j{K) = 
{m + $s} m6 z and the corresponding set of (normalized) cigenfunctions is given 

by 

<t> m (9) := (27T)- 1 / 2 e -*[(m+*B)0-/,fa oo (r)«Jr] _ 



Decomposing into this complete orthonormal set of L 2 ((0, 2ir)), we thus get 
00 CD ( rdr dr ^ r 2 16^ ^ ^ 



Here the symbol ~ stands for the unitary-equivalence relation, the ordinary dif- 
ferential operators on the right hand side should be understood as the Friedrichs 
extensions in L 2 ((0, 00), r dr) of these operators defined initially on Cq°((0, 00)) 
and considered as acting on the respective space L 2 ((0, 00), r dr) ® {4> m } by 
natural isomorphisms. 

It follows that is unitarily equivalent to the harmonic-oscillator Hamil- 
tonian Lho from ([Ml) if $5 £ Z (this can be also understood in terms of a 
gauge invariance). Recall that the form domain of Lho coincides with D{l s ) 
given in (|32|). 

On the other hand, the decomposition (|42t immediately yields the Hardy- 
type inequality 

Loc >y^ with p := dist($ B ,Z) (43) 

in the sense of quadratic forms, which becomes non-trivial whenever $s ^ Z. 
This inequality can be also obtained as a consequence of the Laptev- Weidl 
Hardy- type inequality jTTJ Thm. 3] , representing an elegant modification of © 
for the Aharonov-Bohm field. 

Here we use (|43|) just to conclude with 



®{l o)=W 1 ' 2 (R 2 )nL 2 (M 2 ,\y\ 2 dy)nL 2 (m 2 ,\y\- 2 dy) if $ B ^Z. (44) 

Indeed, if $s ^ Z, the inequality (|43|) enables one to handle the mixed terms 
after developing the second norm in (flU|) and to show that the norm || • is 
actually equivalent to that induced by the form ||Vi>|| 2 + ||yv|| 2 + ||M~ 1 f|| 2 . ^ 
follows that the elements of 53(Zoo) have to vanish at the origin if $b Z (while 
for <E>b S Z they can even diverge there). 

In any case, the form domain S(Zqo) is compactly embedded in L 2 , so that 
the operator has compact resolvent and purely discrete spectrum. The 
latter is explicitly computable as we show now. 

Proposition 3. The spectrum of is given by the discrete set 

, 1 + |m + $ B | 



cr(L c 



nGN, in€2 



The corresponding set of (unnormalized) eigenf unctions is given in polar coor- 
dinates by 

r |m+* B | e - 2 /8£|m+* B |( r 2 /4) ^ _ (45) 

where denotes the generalized Laguerre polynomial. 
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Proof. As a consequence of the decomposition (|42p , we get 



where L™ denotes the ordinary-differential operator on the right hand side 
of (T4"2")) . Consequently, the problem reduces to the study of the spectrum of 
each L™. Fix m G Z. 

The eigenvalue problem L^ip = -E^ leads to an ordinary-differential equa- 
tion that can be cast into Whittaker's equation [TJ 13.1.31]. The latter ad- 
mits solutions in terms of confluent hypergeometric functions M and U, cf [1, 
Sec. 13.1]. Checking the asymptotic behaviours of the functions around the ori- 
gin and at the infinity, it is straightforward to verify that none of the solutions 
belongs to because of singularities, unless the following quantization 

condition for the eigenvalues 

i^ = n + 1 + |W 2 + H „ £ N, {0,1,2,...}, 

holds true. Under this condition, the confluent hypergeometric functions re- 
duce to generalized Laguerre's polynomials [TJ Sec. 22] and the latter lead to a 
complete orthonormal system of eigenfunctions for each . 

More straightforwardly, it follows from fj] 22.6.18] that (f4"5j) is a solution 
of L™ip — E n tp for each fixed m G Z. Then it remains to recall the classi- 
cal fact (see, e.g., Sec. 6]) that {£^}„ £ n is a complete orthogonal set in 
L?((Q,oo),x> i e- x dx). □ 

Remark 3. Note that the spectrum of is doubly degenerate if &b G (^Z)\Z. 
The situation $s G Z corresponds to the classical partial wave decomposition, 
when all the eigenvalues of K are doubly degenerate except for one. Finally, 
the spectrum of Loo is simple in the generic case $s ^Z. 

3.4.2 The strong-resolvent convergence 

Now we are in a position to prove the main auxiliary result of this paper. 

Proposition 4. Let $s g" Z. Then the operator L s converges to in the 
strong-resolvent sense as s —> oo, i.e., 

VFgL 2 , lim WL^F-L^FW = 0. 

Proof. It follows from (|3"6")l that belongs to the resolvent set of L s for all s > 0. 
For any fixed F G L 2 , let us set ip s '■= L^F. In other words, ip s satisfies the 
resolvent equation 

Wg£(Z s ), l,(«,^) = (">F)> (46) 

where Z s (-, •) is the sesquilinear form associated with (I3"2l and (•, •) denotes the 
inner product in L 2 . In particular, choosing %p s for the test function v in 
we have 

^] = W> S ,F)<||^||||F||. (47) 
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Noticing that (j5o]) yields a Poincare-type inequality l a [i/) a ] > IUV'II 2 for anv 
ip G ®(£ s )> we get from (14"T1) the uniform bound 

Hs\\<2\\F\\. (48) 

At the same time, expressing the form l s in polar coordinates (c/ (|T9")) ). the 
bounds (J47]) and (ggj yield 



||a r ^|| 2 < 2 ||F|| 2 , Hr- 1 ^ - *<m/> s )|| 2 < 2 ||F|| 2 , ||W> S || 2 < 32 ||F|| 2 , 

(49) 

where a s (r, 9) :— a(e s ^ 2 r, 9). 

Now, to get a similar estimate on the angular derivative of ip s , we employ 
the Hardy- type inequality © due to Laptev and Weidl as follows. Defining a 
new function u s 6 L 2 by ips(y) = e s ^ 2 u s (e s ^ 2 y) (cf the self-similarity trans- 
formation (O), making the change of variables x = e s l 2 y and recalling we 
have 



l s [ip s }> / V^ s (y) - iA s {y) ips{y)\ dy 
Jr 2 

i 1 2 

Vw s (#) — iA{x) u s (x) dx 



>e*cl ^tdx, 

_ s , I l2 d v- ( 50 

Here the first inequality follows just by neglecting the harmonic-oscillator po- 
tential term in (|32j) . The second inequality is the Hardy inequality ([6]) (the 
hypotheses of [T71 Thm. 1] are clearly satisfied for B £ Cq(R 2 ) assumed in the 
present paper) and holds with a positive c if $>b ^ Z. Using this bound together 
with (J37J) and (gSJ), we therefore get 



As a consequence of (fSTj). we get 

||r- 1 ta.^ a || a < -C\\Ff with C := sup ^-^-a(p,8) 2 . 

c pe(o,oo),ee[o,27r) P 

The finiteness of the constant C follows from the asymptotic behaviour of a(p, 9) 
for small and large p\ indeed, a(p,9) equals the bounded function aoo{9) for all 
sufficiently large p and the bound \a(p,9)\ < /0 2 ||B||oo holds for all p > and 
9 S [0, 27r). It means that we can conclude from the central inequality in (|4"9")l 
the desired uniform boundedness of the angular derivative 

Wr-'dg^W^CWFW, (52) 

where the constant C depends on c and C. 

It follows from (|4"5|) . (|49p and (1521 that {V's}s>o is a bounded family in the 
space £>o := W 1 ' 2 (R 2 ) fl L 2 (M 2 , |2/| 2 rfj/) equipped with the norm || • \\x> given 
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by HV'llso : = II^V'll 2 + Wvi'W 2 + ll^ll 2 - Therefore it is precompact in the weak 
topology of Dq. Let ipoa be a weak limit point, i.e., for an increasing sequence 
of positive numbers {s„}„ e N such that s n — > oo as n — > oo, {?/; Sri }„ e N converges 
weakly to V'oo in ®o- Actually, we may assume that the sequence converges 
strongly in L 2 because £>o is compactly embedded in L 2 . 

Since {ips n }neti converges strongly to ipoa in L 2 and the function e~ s + \y\ 2 
converges locally uniformly for y £ R 2 \ {0} to \y\ 2 as s — > oo, we have 

V*CC^\(0}), ^ 7 =^)_,(^). 

Since C °°(R 2 \{0}) is dense in L 2 (R 2 \{0}) = L 2 (R 2 ) and the uniform bound (J5JJ) 
holds true, we may conclude that 

— / = > i — r m L ■ 

yje- s + | • | 2 | • 

In particular, this means that Voo e L 2 (R 2 ,|yr 2 dy). 

Summing up, we have proved that {^ Sn } n gfj converges weakly to ipoo in the 
space 

Soo :=W 1 - 2 (R 2 )nL 2 (R 2 ,\y\ 2 dy)nL 2 (R 2 ,\y\- 2 dy) 

equipped with the norm || • given by HV'llfo^ := IIV'llli + IIM^Vll 2 - In view 
of (|33]l. coincides with 2)(7oo) as a vector space, so that we know that i/j^ 
belongs to the form domain of the operator L^. Taking v G C^°(R 2 \ {0}) as 
the test function in ([46 j) . with s being replaced by s n , and sending n to infinity, 
we then easily check with help of (|3"5|) that 

ZooKV'oo) = (v,F) . 

Since, by definition (|4"Tj) . C^°(R 2 \ {0}) is a core for 1^, we conclude that ipoo — 
L^F, for any weak limit point of {"0s}s>o- The proof is concluded by recalling 
that we also have the strong convergence of {ips n }«eN to V'oo in L 2 as n — > 00. □ 

Remark 4. The crucial step in the proof is certainly the usage of the Hardy 
inequality in the second inequality of (|50[) . Indeed, it enables one to estimate 
the angular derivative from the central term in (|49[) by controlling the function 
r~ 1 ia s 'ifj s via a uniform bound in s. 

3.4.3 Spectral consequences 

In general, the strong-resolvent convergence of Proposition H] is not sufficient to 
guarantee the convergence of spectra. However, in our case, since the spectra 
are purely discrete (i.e., the operators L s and L M have compact resolvents), the 
eigenprojections converge even in norm (c/ [22] V 

In particular, A(s) converges to the lowest eigenvalue of L M as s — > 00. 
Recalling Proposition [3] and the definition of (3 in (|4"3")l . together with Proposi- 
tion[2]for the case $s € Z, we thus conclude this long subsection by the ultimate 
result: 

Corollary 2. One has 

A(oo) := lim A(s) = (1 + ^/2. 
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3.5 A lower bound to the decay rate 



We come back to (pTlj) . It follows from Corollary [5] that for arbitrarily small 
positive number e there exists a (large) positive time s e such that for all s > s s , 
we have A(s) > A(oo) — e. Hence, fixing e > 0, for all s > s £ , we have 



X(r)dT<- A(r)dr-[A(oo) - e](s - s e ) 
Jo 

< [A(oo) — e]s E — [A(oo) — e]s , 

where the second inequality is due to the fact that A(s) is non-negative for all 
s > (it is in fact greater than or equal to 1/2, cf ([5rj|) ). At the same time, 
assuming e < 1/2, we trivially have 

A(r) dr < < [A(oo) - e]s £ ~[X(oo) - e]s 

also for all s < s e . Summing up, for every s £ [0, oo), we have 

\\u(s)\\ K <C e e-M°°^ s \\u \\ K , (53) 

where C e := e Se > e [H™)-eW. 

Now we return to the original variables {x,t) via (1231) . Using (|2"5"|) together 
with the point-wise estimate 1 < K, and recalling that u = u , it follows 
from d23]) that 

\\u(t)\\ = ||^)|| < ||fi( a )|| x < C £ (1 + 1)-^-] hok 
for every t 6 [0, oo). Consequently, we conclude with 

\\<,-H B t\\ - \\ u {t)\\ ^ r< fi j_ + \-[A(oo)-e] 

ll e |l 2 (k)^l 2 = sup -|| — — <O e (l + t) 
u ei 2 (if)\{o} II u o||k 

for every t G [0, oo). Since e can be made arbitrarily small, this bound implies 

IB > A(oo) = (l + /3)/2. 

This together with Corollary [T] proves Theorem [TJ 



3.6 A global upper bound to the heat semigroup 

Theorem [T] provides quite precise information about the extra polynomial decay 
of solutions u of ([T]) in the sense that the decay rate jb is better by a factor /3/2 
with respect to the field-free case. On the other hand, we have no control over 
the constant C 7 in (J3J (in principle it may blow up as 7 — > 7s). We therefore 
conclude this section by establishing a global (in time) upper bound to the heat 
semigroup (i.e. we get rid of the constant C 7 ) but the prize we pay is just a 
qualitative knowledge about the decay rate. It is a consequence of (|3"Tj) and the 
following improvement upon f|36[) : 

Proposition 5. Let B ^ 0. Then 

Vse[0,oo), A(s)>l/2. 
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Proof. The fact that A(s) > 1/2 for all s > 0, regardless of the presence of B, is 
due to (1551) . To show that the inequality is strict, let us assume by contradiction 
that B ^ and that A(s) = 1/2 for some s > 0. Let w denote a corresponding 
eigenfunction of L s . We have 

o = i«M - \ \H 2 > p\v\\\ 2 + i ||HI 2 - 5 IN 2 > 0, (54) 

where the first estimate is due the diamagnetic inequality (|35[) and the second 
inequality follows from the variational characterization of the lowest eigenvalue 
of Lho (which is 1/2). Hence, \v\ must coincide with the ground state of £ho- 
The latter is unique (up to a normalization factor) and can be chosen positive. 
Consequently, v as an eigenfunction of L s is unique (up to a normalization 
factor) and can be chosen positive. The equality in (|54l) then implies 

\\A s v\\ =0. 

Hence, employing the continuity of A s and v, A s = identically, which in turn 
implies that B is trivial, too. □ 

Remark 5. If B ^ but $s 6 Z, then it follows from Proposition [2] that A(s) 
converges to 1/2 as s — » oo, despite of the strict inequality of Proposition [5] 

Combining Proposition [5] with Corollary^ we see that the number 

c B := inf A(s)-l/2 (55) 

se[o,oo) 

is positive if, and only if, ^ Z. In any case, (1511) implies 

||«(»)||iC < \\u \\Ke-^ +1 / 2 > 

for every s £ [0, oo). Using this estimate instead of (j53|) . but following the same 
type of arguments as in Section 13.51 below ((53]) , we thus conclude with: 

Theorem 2. Let B G C^(R 2 ). We /iawe 

Vf G [0, oo), ||e- ffB *|| L2(i ^ L2 < (1 + t)-< c « +1 /2) , 
w/iere c s > i/ $s ^ Z /and c_b = i/ $s G Zj. 

4 Conclusions 

The results of this paper show that the presence of a local magnetic field leads 
to an improvement of the large-time decay of solutions of the heat equation. As 
explained in the introduction, the physical reason behind the better decay is the 
diamagnetic effect of the magnetic field. This is mathematically expressed in 
terms of the existence of Hardy-type inequalities for the magnetic Schrodinger 
operator, which played a central role in our proofs. 

It follows from our analysis that the solutions of the heat equation for large 
time behave as if the magnetic field were replaced by the singular Aharonov- 
Bohm magnetic field with the same total flux. This also explains the role of the 
"degeneracies" $b 6 Z, since they correspond to the flux quanta for which the 
Aharonov-Bohm effect is not detectable in our setting (c/ [TSJ Sec. 10.2.4]). 
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Since no assumptions about the symmetry of the magnetic field are made 
whatsoever in the present paper, it can be also viewed as a generalization of 
some of the results of the recent work [T3] 

It is known that the method of self-similarity variables yield sharp decay 
rates. Therefore, we conjecture that the inequality of Theorem[T]can be replaced 
by equality, i.e., jb — (1 + j9)/2 for any field B. A careful analysis in the spirit 
of [9J [7] could even show that the solutions can indeed be approximated by the 
solutions of the effective equation with the Aharonov-Bohm field. 

The question of optimal value of the constant cb (and its quantitative de- 
pendence on the magnetic field B, especially on the total magnetic flux 
from Theorem [2] also constitutes an interesting open problem. Note that the 
constant is related to the decay rate bycB + l/2<7 B . 

We expect the same decay rates if the assumption about the compact support 
of B is replaced by its fast decay at infinity. However, it is quite possible that 
a slow decay of the field at infinity will improve the decay of the solutions even 
further. In particular, can jb be strictly greater than (1 + (3)/2 if B decays to 
zero very slowly at infinity? 

The characteristic hypothesis of the present paper, under which the addi- 
tional decay rate is proved, is that the total magnetic flux is not an integer. Even 
if this hypothesis is violated and the magnetic field is non-trivial, however, there 
exist Hardy- type inequalities for the magnetic Hamiltonian |21j . Moreover, it 
is shown in [T3] that there is a logarithmic improvement to the decay of the 
heat kernel even if the total flux is an integer (and the field is rotationally sym- 
metric) . It should be possible to establish the additional logarithmic decay rate 
by the method of the present paper (without any symmetry assumption), by 
simply computing the next term in the asymptotic expansion of the eigenvalue 
A(s) as s — > oo. 

More generally, recall that we expect that there is always an improvement 
of the decay rate for the heat semigroup of an operator satisfying a Hardy-type 
inequality (cf [HI Conjecture in Sec. 6] and Q21 Conjecture 1]). The present 
paper confirms the general conjecture in the particular case of two-dimensional 
magnetic Schrodinger operators. 

Finally, it would be interesting to examine the effect of the presence of the 
local two-dimensional magnetic field in other physical models. As one possible 
direction of this research, let us mention the problem of Strichartz estimates for 
the magnetic Schrodinger equation studied recently in higher dimensions in [B]. 

Acknowledgement 

The author would like to thank Carlos Mora-Corral for useful discussions. The 
work was partially supported by the Czech Ministry of Education, Youth and 
Sports within the project LC06002 and by the GACR grant P203/11/0701. 

References 

[1] M. S. Abramowitz and I. A. Stegun, eds., Handbook of mathematical functions, 
Dover, New York, 1965. 

[2] G. E. Andrews, R. Askey, and R. Roy, Special functions, Encyclopedia of Math- 
ematics and Its Applications, vol. 71, Cambridge University Press, 1999. 



19 



[3] A. Balinsky, A. Laptev, and A. Sobolev, Generalized Hardy inequality for the 
magnetic Dirichlet forms, J. Stat. Phys. 116 (2004), 507-521. 

[4] H. Brezis, Analyse fonctionnelle: Theorie et applications, Dunod, 2002. 

[5] X. Cabre and Y. Martel, Existence versus explosion instantanee pour des 
equations de la chaleur lineaires avec potentiel singulier, C. R. Acad. Sci. Paris 
329 (1999), 973-978. 

[6] P. D'Ancona, L. Fanelli, L Vega, and N. Visciglia, Endpomt Strichartz estimates 
for the magnetic Schrddmger equation, J. Funct. Anal. 258 (2010), 3227-3240. 

[7] G. Duro and E. Zuazua, Large time behavior for convection-diffusion equations 
in M. N with asymptotically constant diffusion, Commun. in Partial Differential 
Equations 24 (1999), 1283-1340. 

[8] M. Escobedo and O. Kavian, Variational problems related to self-similar solutions 
of the heat equation, Nonlinear Anal.-Theor. 11 (1987), 1103-1133. 

[9] M. Escobedo and E. Zuazua, Large time behavior for convection-diffusion equa- 
tions in R N , J. Funct. Anal. 100 (1991), 119-161. 

[10] W. D. Evans and R. T. Lewis, On the Rellich inequality with magnetic potentials, 
Math. Z. 251 (2005), 267-284. 

[11] S. Fournais and B. Helffer, Diamagnetism, Spectral Methods in Surface Supercon- 
ductivity, Progress in Nonlinear Differential Equations and Their Applications, 
vol. 77, Birkhauser Boston, 2009, pp. 19-30. 

[12] M. Frass, D. Krejciffk, and Y. Pinchover, On some strong ratio limit theorems 
for heat kernels, Discrete Contin. Dynam. Systems A 28 (2010), 495-509. 

[13] D. J. Griffiths, Introduction to quantum mechanics, Prentice Hall, Upper Saddle 
River, NJ, 1995. 

[14] H. Kovaffk, Heat kernels of two dimensional magnetic Schrddinger and Pauli 
operators, preprint on http://calvino.polito.it/~kovarik/ 

[15] D. Krejciffk and E. Zuazua, The asymptotic behaviour of the heat equation in 
a twisted Dirichlet- Neumann waveguide, J. Differential Equations, to appear; 
doi:10.1016/j.jde. 2010. 11.005. 

[16] , The Hardy inequality and the heat equation in twisted tubes, J. Math. 

Pures Appl. 94 (2010), 277-303. 

[17] A. Laptev and T. Weidl, Hardy inequalities for magnetic Dirichlet forms, Oper. 
Theory Adv. Appl. 108 (1999), 299-305. 

[18] K. Pankrashkin and S. Richard, Spectral and scattering theory for the Aharonov- 
Bohm operators, Rev. Math. Phys., to appear; preprint on arXiv:0911.4715 [math- 
ph]. 

[19] M. Reed and B. Simon, Methods of modern mathematical physics, IV. Analysis 
of operators, Academic Press, New York, 1978. 

[20] J. L. Vazquez and E. Zuazua, The Hardy inequality and the asymptotic behaviour 
of the heat equation with an inverse-square potential, J. Funct. Anal. 173 (2000), 
103-153. 

[21] T. Weidl, A remark on Hardy type inequalities for critical Schrddinger operators 
with magnetic fields, Oper. Theory Adv. Appl. 110 (1999), 345-352. 

[22] J. Weidmann, Continuity of the eigenvalues of self-adjoint operators with respect 
to the strong operator topology, Integral Equations and Operator Theory 3 (1980), 
no. 1, 138-142. 



20 



